The International Trade Network: weighted network analysis and modelling 
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Tools of the theory of critical phenomena, namely the scaling analysis and universality, are argued 
to be applicable to large complex web-like network structures. Using a detailed analysis of the real 
data of the International Trade Network we argue that the scaled link weight distribution has 
an approximate log-normal distribution which remains robust over a period of 53 years. Another 
universal feature is observed in the power-law growth of the trade strength with gross domestic 
product, the exponent being similar for all countries. Using the 'rich-club' coefficient measure of 
the weighted networks it has been shown that the size of the rich-club controlling half of the world's 
trade is actually shrinking. While the gravity law is known to describe well the social interactions 
in the static networks of population migration, international trade, etc, here for the first time we 
studied a non-conservative dynamical model based on the gravity law which excellently reproduced 
many empirical features of the ITN. 
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Are the large real or man-made networks with com- 
phcated and heterogeneous connections scale invariant 
and obey the well known scaling analysis and universality 
concepts of Statistical Physics? In the recent years, ex- 
tensive research effort have been devoted to analyzing the 
structure, function and dynamics of complex networks 
relevant to multi-disciplinary fields of science [H, 0, 
Indeed the scale-free networks reflect such scale invari- 
ance in the link-node structures of the electronic com- 
munication networks like the Internet, world-wide web, 
protein interaction networks and even in research collab- 
oration networks. More recently it has been observed 
that strengths of links of the networks, as called weights 
in graph theory, also have very interesing properties and 
can shed much light into the understanding of the details 
of the network [J 0, Q • Lately it has been proposed 
that the International Trade Network (ITN) i.e., the sys- 
tem of mutual trading between different countries in the 
world can also be viewed as an interesting example of 
real- world network 0, [O, El Q lM] • 

In this paper we study the ITN as an excellent ex- 
ample of the weighted networks obeying the scale invari- 
ance and universality properties where the extent of trade 
between a pair of countries can be treated as the link 
weight p^. Il7||. It is observed that the suitably scaled 
link weight distribution over many years can be approx- 
imated well by a log-normal law. The nodal strength 
measuring the total volume of trade of a country is seen 
to depend non-linearly on the country's Gross Domestic 
Product with a robust exponent. Many features observed 
in the analysis can be explained with a simple dynam- 
ical model, which has the well-known grav ity model of 
international trade as its starting point [18]. 

In the ITN, a node depicts a country and an undi- 
rected link exists between any pair of nodes if the trade 
volume between the corresponding countries is non-zero. 
Both the number of nodes N and the number of links 
L in the ITN's show annual variation and grow almost 

systematically over 53 years, e.g. in 1948 N — 7Q and 
T — ^AQA „7V,£ir-Qoc! ir. onnn at — 187 cnri T — ^^\0K0 



On the other hand the link density L/ [N{N - l)/2] is 
observed to remain roughly constant with values around 
0.52 for the same period. The annual trade data are ex- 
pressed in millions of dollars (Af $) of imports and exports 
between countries i and j using four different quantities 
expjj , expj j, impy and imp^ j [H, [i3]- -Due to differences 
in the reporting procedures, there are usually small de- 
viations between exports exp^^ from i to j and imports 
impjj to j from i. Therefore, we define the link weight 
as a measure of the total trade volume between the 



two countries (in M$), as follows 



Wi 



(exp, + exp j + imp, + imp ,)/2. 



(1) 



This quantity tends to average out the aforementioned 
discrepancies. The distribution of weights is observed to 
be broad with the smallest non-zero trade volume being 
less than 1 M$ and the largest of the order of 10^A/$. 
The number of links with very small weights is quite 
large, whereas only a few links with very large weights 
exist (Fig. 1). The tail of the distribution consists of links 
connecting very few high-income countries [l9j . The av- 
erage weight per link is observed to grow during the in- 
vestigated period, from 15.54 M% in 1948 to 308.8 M$ in 
2000. 

The probability density distribution of the link 
weights is defined as the probability PToh{w)dw that a 
randomly selected link has the weight between w and 
w + dw. The probability decays systematically with in- 
creasing weight and has a long tail with considerable fluc- 
tuation. We note that inferring the form of a broad prob- 
ability distribution based on a relatively small number of 
data points has its difficulties; in the case of ITN, log-log 
plots of the distributions (not shown) display small inter- 
mediate linear regions which are (too) often interpreted 
as power laws. However, the tails of the distributions 
show clear curvature and much better fits to data are 
obtained by using a log-normal distribution of the form 
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FIG. 1: (Color online) A subnetwork of the ITN for the year 
2000, where only links with the highest 4 % of weights and 
the associated nodes (countries) are included [lO] , yielding in 
total 80 nodes and 411 links. The node size is proportional 
to its strength and hnk color to its weight. Link weights are 
defined as the volume of annual trade between two countries 
in M$. 

where the constants defined as wq = exp((ln(ti;))) and 
a = {((ln(i(;))^) — (In(w))'^}^/^ are observed to have differ- 
ent values for the weight distributions for different years. 
However, one can get a plot ind epend ent of wq and a 
by drawing —2a^ ln[Prob{ln(K;)}A/27ro^] as a function of 
ln(ti;/i/7o), which gives a simple parabola y = for all 
years (Note that Prob{ln(it;)}(i{ln(w)} — 'Pmh{w)dw im- 
plies Prob{ln(i(;)} = wProb(w)). In Fig. 2 we show such 
a plot, where the data has been aggregated for the five- 
year periods to reduce noise, i.e., 1951-55, 1956-60, 
1996-2000, each represented with its own symbol. It is 
clearly seen that the data points for each period fall close 
to the parabola y = x"^, displayed as the solid line. There- 
fore, we conclude that the annual weight distributions 
are reasonably well approximated by the log-normal dis- 
tribution. This log-normal distribution of trade volumes 
has been discussed also in [l^ . Note also that the trade 
imbalances have been reported to follow a log- normal dis- 
tribution [Tsj. 

For weighted complex networks, the strength Si of 
node i is defined IJI, which in the case of 

ITN corresponds to the total volume of annual trade as- 
sociated with a node. Intuitively, one can expect that in 
general the strengths of high- GDP. countries are higher 
than those of the low-GDP countries. To see this in a 
quantitative fashion, we have utilized an elastic constant 
7 to measure how changes in strength respond to changes 
in GDP (we have used total Gross Domestic Product of 
the country and not the GDP per capita without tak- 
ing into account of inflation). Hence, we define 7 as: 
[ds / s) / [dG / G) = 7, which results in a power-law rela- 
tionship Si cx Gj . In Fig. 3 we plot the strengths Si vs. 
GDP Gi for 22 different rmintrie.'j- renresentinp' a, prnnd 
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FIG. 2: (Color online) The probability distribu- 
tion of link weight s of the ITN. This plot displays 
-2cr^ ln[Prob{ln(?j;)}V27rcr2] as a function of \n{w / wa) . 
The data points scatter around the parabola y — (solid 
line) , which implies a log-normal probability distribution (see 
text). Different symbols depict empirical distributions where 
the data have been aggregated over five-year periods. The 
inset shows the distribution of all link weights over the time 
span of 50 years from 1951 to 2000. 



mix of economic strengths, i.e., 10 high, 5 middle and 7 
low income countries. It is observed that the strength of 
each country grows non-linearly with its GDP with ap- 
proximately the same slope. In the inset of this figure 
we show the probability distribution of 7 values for 168 
different countries. The distribution has a long tail and 
the 7 values of 12 countries are found to be larger than 
2. A detailed inspection reveals that majority of these 
12 countries are those originated after Soviet Union, Yu- 
goslavia and Czechoslovakia were fragmented. An overall 
average of the growth exponents has been estimated to be 
7 = 1.26 which comes down to 1.06 if these 12 countries 
are not considered for the averaging. The peak value of 
the distribution occurs very close to 7 = 1. Interestingly, 
Irwin [2l[ has observed earlier that the total world ex- 
port volume varies as a function of the total world real 
GDP to the power of 1.16, along with other factors. In 
comparison, however, our observation reveals a more de- 
tailed picture, indicating that the total trade volumes of 
individual countries are also approximately power laws as 
function of their individual GDPs, with exponents close 
to this value. 

Non-trivial correlations among the nodes of differ- 
ent real- world networks have been observed. For an 
unweighted network this means that large degree nodes 
are connected among themselves forming a club. More 
precisely such a club consists of a subset of Uk nodes 
whose degree values are at least k. A rich-club coefficient 
(RCC) is measured as (p{k) = 2Ek/[nk{nk — 1)] where 
Ek is the number of links actually exists in the club and 
[nfc(nfc — 1)]/2 is the total number of node pairs in the club 
j22j . A high value of 4>{k) implies that members are in- 
deed tightly connected. However it has been realized that 
nnlv this definition i.s not enmicrh. .since with this measure 
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FIG. 3: (Color online) The dependence of strength Si{Y) of 
country i measuring its total volume of trade on its GDP 
Gi{Y) in year Y, plotted for each country (color coded) over 
the 53 year span from 1948 to 2000. General trend corre- 
sponds to a non-linear growth with an average exponent of 
around 1.2. The inset shows the distribution of exponents 
associated with individual countries. 



even uncorrelatcd random graphs show some rich-club ef- 
fect as well [2^]. It is suggested that one needs to define 
a 'null model' or the maximally random network (MRN) 
keeping the nodal degree values {ki} preserved, measure 
the corresponding RCC 4>ran{k) and observe the variation 
of ratio p{k) = (j){k) / (firanik) . We have executed the same 
analysis for the ITN of the year 2000 and generated the 
MRN using the pairwise linkend exchange method [2^ . 
However it is observed that the variations of (/)(fc) and 
4'ran{k) with k are nearly same and consequently p{k) is 
nearly equal to unity for the whole range of degree values. 

Next we studied the rich-club effect of the same ITN 
but now considering it as a weighted network. The rich- 
club is now defined as the subset of nodes whose strengths 
are at least s controlling a major share of the world's 
trade dynamics. The RCC of the weighted network is 
defined as: 



Rw{s) = 2Y,i^jWij/[ns{ns - 1)], 



(3) 



The corresponding maximally random weighted network 
(MRWN) has been generated keeping both the nodal de- 
grees {ki} as well as the nodal strength values {si} pre- 
served. 

To generate the MRWN from the MRN of the ITN of 
the year 2000 described above a self-consistent iteration 
procedure is used to obtain the link weight distribution 
consistent with the nodal strength list {si}. We start 
assigning arbitrary random numbers as the weights Wij 
to all links maintaining that the weight matrix is always 
symmetric i.e., Wij — Wji. For an arbitrary node i the 



difference 6i 



is calculated. Weights of all hi 



links meeting at the node i are then updated as 

Wij Wij + Si {wij/Y^jWij) . 



(4) 




to balance Si and T.jWij. By repeated iterations the link 
weifrhts niiicklv mnvercenpe and attain mnsistencv with 



FIG. 4; (Color online) (a) The weighted rich-club coefficients 
R^{s) (black) and ii™"(s) (red) for the ITN of the year 2000, 
(b) The ratio pran{s) of these two coefficients is nearly equal 
to unity over the whole range of variation and (c) fm{s), of 
the international trade volume among the rich-club members 
having strengths s and above. 



nodal strengths {si}. It is checked that (siSj) ^ Wij 
relation is well satisfied for this MRWN. 

In Fig. 4(a) we plot both i?^(s) for the ITN and 
i?™"(s) for the MRWN with the scaled nodal strength 
s/smax- The two measures are found to be nearly same, 
grow like gO.ss ^g^j-gg g values and their ratio Pran{s) 
is nearly equal to one except for a few values of s near 
S^ax (Fig. 4(b)). 

To explain this we observe that only 15% elements 
of the adjacency matrices of the ITN and the MRN are 
different. Therefore a typical node of ITN retains the 
links to most of its neighbors even after maximal ran- 
domization. This is because of the high value of the link 
density (59%) of the ITN for the year 2000. As a result 
Pran{k) as wcU as Pran{s) are nearly equal to unity. This 
implies that the pairwise link connections and the asso- 
ciated link weights of the original ITN are very close to 
those of the corresponding MRWN. In fact one can say 
that the original ITN is a typical member of the different 
random configurations of the MRWN when the {fc^} as 
well as {si} sets of the ITN are preserved. These results 
on the randomized unweighted ensemble confirms the re- 
sults in (lo| when compared with randomized networks 
with fixed degree sequence as in [14 1 . 

Zhou and Mondragon [23| observed a very similar 
behavior for the Internet statistics. Following them we 
conclude that rich nodes in the original ITN and in its 
cnrresnnndinp' MRN and MRWN are ticrhtlv cnnnecteH 



4 





FIG. 5: Model results: (a) Scaled weight distribution fits 
well with the simple parabola y = . (b) The GDPs of 
the individual countries have a broad distribution, which in 
the tail region seems to follow roughly a power law behavior 
(shown as a line with exponent ~1.92). (c) The strength 
correlation (siSj) as a function of link weight is shown to 
grow approximately as (stSj) oc Wif^. 



and the similarity of the rich-club connectivity in the ITN 
structure (with and without weights) does not imply that 
ITN lacks a rich-club structure. 

In fact the presence of the rich-club effect is evident 
even if we simply analyze the variation of the fraction 
/to(s) of the total volume of trade taking place among 
the members of the club. This is depicted in Fig. 4(c) 
which shows that fw{s) remains close to unity until high 
values of s/smax, then decreases gradually to 1/2 at 
s/smax ~ 0.11, and finally drops sharply. Very interest- 
ingly we also observe that the size of the rich club con- 
taining 50% of the total volume of annual trade shrinks 
almost systematically from 19% in 1948 to 8% in 2000. 

The very heterogeneous distribution of trade vol- 
umes in the ITN is also reflected in the average pair corre- 
lation function (siSj) of nodal strengths and in its power 
law dependence on the link weights Wij . Links with high 
weights Wij Wmax obviously must connect pairs of 
nodes of high strength, and for them (siSj) ^ ^max- 
On the other hand for links of weights around unity, 
(siSj) ~ Smax- Assuming that Wmax itself is of the order 
of Smax, we find an upper bound for the exponent v — 1 
describing the variation of {siSj) ~ w^y Our analysis of 
the ITN yields, however, a somewhat smaller value of u 
being between 0.65 and 0.90 for different financial years 

bet^ween 1948 and 9(10(1 Thf> dfinendpnrp nf thp wpicrVit 



distribution on the underlying topological structure of 
the ITN is studied by measuring the average strength of 
a node as a function of its degree, which turn out to ex- 
hibit strong degree of non-linearity: {s{k)) oc /c^ where n 
varies between 3.4 and 3.7 for the same period. 

We now develop a dynamical model based on the 
well-established gravity model [18| used in social and 
economic sciences to describe the flow of social interac- 
tion between two economic centers i and j as a func- 
tion of their economic sizes , rrij and distance £ij : 
Fij = Gminij/ifj. This equation has been generalized 
to the parametric form [24] 

^^.-<(^/Sm»^)- (5) 

where the exponents a and (3 usually range between 0.7 
and 1.1 where as 6 is observed to be around 0.6 [23 |. 

In our model, we assume a unit square to repre- 
sent the world and N points distributed at random po- 
sitions representing the capital cities of different coun- 
tries. Initially the GDP values mi{i — l,N) are ran- 
domly assigned with uniform probability such that the 
total GDP is unity: Tii^i^^rrii = 1. Then we let the dy- 
namics start, which is essentially a series of pairwise in- 
teractions. There are different pairwise wealth exchange 
models studied in the literature which calculate the dis- 
tribution of wealths in a society [1^ . A pair of countries 
is randomly selected (1 < i,j < N) for a trans- 
action and time t is measured as the number of such 
transactions. In a transaction the selected countries in- 
vest the amounts Fij and Fji calculated using Eq. 
Then the total amount of investment Fij = Fij + Fji is 
randomly shared between the two countries as a result of 
this transaction, as follows: 

m, = m,j - Fij + ePij + A,, (6) 
ruj = ruj ~^Fji + {l- e)Fij + Aj . (7) 

Here e is a random fraction freshly drawn for every trans- 
action. The random sharing of Fij is justified by the fact 
that while the Gravity law describes the average inter- 
action in terms of the strengths and distances of separa- 
tion, the actual amount of trade depends on other fac- 
tors, many of them are political. In this idealized model 
countries are not allowed to make debt, which in turn 
makes the dynamics non-conservative through the pa- 
rameters Ai and Aj. It holds for these parameters that 
Ai = a Fij < rrii and Aj = if Fji < rrij. How- 
ever, if for some transaction Fij > m.^ or Fji > rrij then 
we add = Fij — rrii or Aj = Fji — rrij such that af- 
ter the transaction, GDP balance does not become neg- 
ative. Also after the transaction the individual GDP's 
are rescaled rrii rrii/'Sjmj for the total GDP to re- 
main unity. It is observed that a large number of pair- 
wise transactions leads to a stationary state where (m^) 
fluctuates with time around a steady mean value. Any 
time after reaching the stationary state, the dynamics is 
used to construct a model ITN such that links are es- 
tablished between c.nimtrips i and i whpnevpr thprp i.s a. 
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transaction between them. We let the dynamics run un- 
til a pre-assigned link density (typically 0.3-0.5) has been 
reached. For example to generate a network correspond- 
ing to the ITN of the year 2000, we take N = 187 and 
continue the exchange dynamics till L = 10252 distinct 
links are dropped corresponding to the link density 0.59. 
The weight of a link is then defined as the total amount of 
investment between pairs of countries in all transactions. 

For comparison the weight distribution of our model 
networks is analyzed in the same way as the real ITN 
data, and it turns out that an excellent consistency with 
the simple parabola y = is observed for parameter val- 
ues a — 1/2,(3 = 1,9 = 1/2, within a tolerance of 0.2 for 
all exponents, see Fig. 5. We also find that the GDP dis- 
tribution is broad showing a power law like behavior for a 
short interval of m (with exponent ^^1.9) before finite size 
effects set in. This is to be compared with the real-world 
data where the GDP distribution of different countries 
has been argued to be consistent with the Pareto law 
[26| . Finally, as in the real ITN, the two-point strength 
correlation is seen to grow like {siSj) oc w'/j with v « 0.98 
in the large weight limit, compared to the range of val- 
ues 0.65 - 0.90 in the real ITN. In addition we find that 
the cumulative degree distribution is consistent with the 



results on the real ITN. 

To summarize, we have shown that the weighted In- 
ternational Trade Network can be looked upon as an 
excellent example of a complex network obeying scale- 
invariance and universality features. The scaled distri- 
butions of annual world trade volumes between countries 
collapse well to a log-normal distribution and it remains 
unchanged over a span of 53 years implying robustness 
or universality. Secondly, the nodal strength measuring 
the total trade volume associated with a country grows 
non-linearly with its GDP with a robust exponent. Also 
a large fraction of the global trade is controlled by a 
club of few rich countries which shrinks its size as time 
goes on. Finally, the main features of the real-world ITN 
have been reproduced by using a simple non-conservative 
dynamical model starting from the well-known gravity 
model of social and economic sciences. 
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